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ABSTRACT 


The main aim of present study is to examine a Non-Newtonian mathematical model of two phased blood flow in human 
pulmonary capillary, keeping in view the nature of pulmonary blood circulation. Herschel Bulkley Non — Newtonian model 
in Bio — fluid mechanical setup is applied with respect to the help of clinical data in case of Asthma for hemoglobin versus 
blood pressure. In present study overall presentation is in tensorial form and the solution technique adopted is analytical 


as well as numerical. 
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INTRODUCTION 


1. Structure and Function of the Lungs 


The human lungs are paired organs in the chest and divided in two lobes. Where the right lung has three lobes, the left one 


has two lobes. It is smaller to accommodate to the heart, which takes up the space of the third lobe in the chest cavity. 


Figure 1 
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Figure 2. 


Figure 3. 


2. Mathematical Modeling (Choice of Frame of Reference) 


Manoj Kr. Srivastava 


A frame of reference on mathematical modeling for moving blood we observed by the three - dimensional orthogonal 


curvilinear co-ordinate system, briefly prescribed as E3,called as 3-dim Euclidean space (tensorial form). 


2.1. Equation of Continuity 


Let the volume portion covered by blood cells in unit volume be X, this X is replaced by H/100, where H is the 


Hematocrit the volume percentage of blood cells. Then the volume portion covered by the plasma will be 1-X. If the mass 


ratio of blood cells to plasma is r then clearly. 
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pe Xpc 


where, Pc and Pp are densities of blood cells and blood plasma respectively. 


The both phase of blood, move with the common velocity. Hence, the equation of continuity for two phase by the 


principle of conservation of mass is 


HP) » xpey),i= 0 
Ot (2) 
And 
OUST? ee Bye 
Ot (3) 


Where, v is the common velocity of two phase blood cells and plasma. Again (Xpcv’) is co-variant derivative of 


(Xpev') with respect to X'. In the same way (1-X)Pc Vv’), i is co-variant derivative of ((1-X)pcV') w.r. to X! 


If we define the uniform density of the blood Pm then we consider 
l+r +r 1 
— + 
Pm = Pe —-~Pp 


(4) 
From equations (2) &(3) we consider 
Op 
m . 
Or + (p mv" di = 0 
t 
(5) 
2.2. Equation of Motion 
The equation of motion for the two phase of blood cells is 
Ovi ( iyi ij jk qi, + 
aa +\XpV/ WV; = — Xp, jg” + Xnc(g” V,), J 
’ (6) 
Similarly, the equation of motion for plasma is 
avi Ss ; ; ae 
(x )Pp— + (0-x)Ppv'W) =-(1-X)p. jg? + 0-X)p. ia? + 0—X)ne VD. j o 


Now, adding equations (6) & (7) and by relation (4), the equation of motion for blood flow with the both phases 
will be 
i . ik ot : 
=~ Pri es Vp 


(8) 
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Where nti ae = X)p is the viscosity coefficient of blood as a mixture of two phase. As the 
velocity of blood flow decreases, the viscosity of blood increases. In this situation, the blood cells line up on the axis to 
build up rouleaux. Hence a yield stress is produced. Though this yield stress is very small, even then the viscosity of blood 
is increased nearly ten times. The Herschel Bulkley law holds good on the two phase blood flow through veins Lungs, 


veinules and whose constitutive equation is 


T=n,e¢° +T, (T'2T,) and 


e=0(T'<T,) where, (T'=> qT, ) 


T, is the yield stress. When strain rate e = 0 a core region is formed 


which flows just like a plug. Let the radius of the plug be r,. The stress acting on the surface of plug will be T,. Equating 


the forces acting on the plug, we get, 
Pa, =T,2 
TU p — 4p TIL, 


= 7 = gee 
ig (9) 


Blood Vessel Ping region or rouleaux 


Fig (1): Herschel Bulkley blood flow 


Figure 4 
The Constitutive equation for test part of the blood vessel is 
T' =nme"+T, OrT’- T, =n,,e° =T, 
Where, Te = effective Stress 
Whose generalized form will be 
T! = — Pg! +T° where TT =n,,(e!)" While e? = g* V. 


Now we describe the basic equation for Herschel Bulkley blood flow then 
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(i). Equation of Continuity 


1/ g( gvi,=0. 


(ii) Equation of Motion 


avi - : 

Dy — + PmvV'V'., =-T, } 

m at e,] 
(10) 


Where, all the symbols have their usual meanings. 


2.3. Analysis 


Since, the blood vessels are cylindrical; the above governing equations have to be transformed into cylindrical co- 


ordinates. 
1 2 3 
xX =r, XX" =0,X =Z, 
Matrix of metric tensor in cylindrical co-ordinates is 


10 0 
Lg] = |0 ro 
00 1 


While matrix of conjugate metric tensor is 


[g*] = 0 


0 
; 

00 1 

Where, as the Chritoffel's symbols of 2™ kind are 


ed Bas 


Remaining others are zero. 


Relation between contra variant & physical components of velocity of blood flow will be 


| 1 = _ wl 
Si a Vv, = v.= Vv 
2 2 
y82v = Vo => Ve=rv 
| 3 3 
833 = Vz =e v= v 
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. . ij oi Page” 
Again the physical components of -p,jg” are - oe 


Equation (9) and (10) are transformed into cylindrical form so as to solve them as power law model then 
u 
dv _| Pr |" 
dr \ Qn, 
dp 


Where, pressure gradient dz 


dy | (ely 


Nn 


= P 


2 2 
dr ue 
by equation (9) 
1 Pa 
av |> peat: 
dr Nin 
(1) 
Substituting the value of Tp from (7) into (11), then 
is 
1 1 a 
dv | 9 pr- Brie 
dr nm 
dv 7 1 
— = —(p/2m)" (r=r, 
dr (12) 
Integrating under the no slip boundary condition: v = 0 at r= R then we have 
Ion : ! 
V=(P/2n,,)r (er, (rr 
This is the formula for velocity of blood flow in Lungs veinules and veins. 
Putting r =r, then the velocity V, of plug flow is 
1 1 
V, = —~ (P2nm)r(R-4, jn" 
n+l (14) 
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2.4. Result 

Observations 

Hematocrit vs. blood pressure from an authorized. Homeo Chikitsa Kendra, Allahabad by Dr. K. K. Mishra. 
Patient: Pyare lal 


Diagnosis: Asthma. 


B-H Table 
Table 1 
S.N Date HB(Hemoglobin) | B.P. (Blood Pressure) | Hematocrit | Hemoglobin 
1 10/9/20 12.2 140/100 36 12 
2 16/9/20 11.6 130/90 35.4 11.8 
3 20/8/20 11.4 130.3/80 33.6 11.2 


The flow flux of two phased blood flow in Lungs, veinules and veins is 
Q= [ 2arv, dr + q 2arvdr 


[rn 2 (2) an a+ 


2nm | 


(me * (rom | ea} Ja 


n+l 
Solving (12) & (14) then 


1 


2mn P "(r +k r a nn P “(r a r(r—r, % (c—r, Jr" 
(n+1)\ 2nm 7 2|, @+tl)\2nm) | 2 ’ 1 (: iE 


Q = 425.136 ml. / min 


R= 1,1)=.999 
Le (Pascal - sec.) 
Nn = 0.027 (Pascal- sec.) 
.H = 24 
H 24 
=n.xX + 1—X) where, X =—— = —=.24 
fj, = Hx Sq, GX) io Te 
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26n* + 33n* + 9n 
6n> + 11ln* +6n +1 


1 
1827.23 = (1728.40)" 


Solved by Numerical method, then n = 1.127 


‘ 26n? + 33n?+9n ) 
AP=(003213H +.0039) (82723) x | ~="" #7" 
6n° +11n~ + 6n +1 


36.1 | 354 | 336 | 33.2 
| 193.312 | 190.193 | 180.843 | 177.721 | 


H 
| AP 


Blood Pressure 


33 33.6 35.4 36 
Hematocrit 


Figure 5 


CONCLUSION 


It is clear that when Hematocrit is increased the Blood pressure also increased. Hence Hematocrit is proportional 


to blood pressure. 
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